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I. Crystal Model 

Consider a cubic crystal formed by a periodic potential V(r) = U(x) + U(y) + U(z) , where 

U(x + a) = U(x) and a is the lattice constant. The corresponding Hamiltonian eigenfunctions 
are the Bloch functions [1] 

^Mk(r) = e*V00 = ^ (x)0 nk> (y)0 qK (z) , (1) 

V ( r ) = u mks (x)u nK (y)u qk (z) ; u mk (x + a) = u mk (x) , (2) 

where ^(r) and u mk (x) are three- and one-dimensional Bloch amplitudes correspondingly, 
M = (m,n,q) is the band index, and k = (k x ,k y ,k z ) is the quasi-wavevector of the state. The 
one-dimensional Bloch functions <j> mk (x) = e' kx u mk (x) are normalized as follows [2] 



GO 

j u * nk (*K* ( x ) dx = -?- S mn and j" f nq (x)(/> mk (x)dx = S mn S(q - k) . 



(3) 



The set of the Bloch functions is complete. The energy of the state (1) is given by 

E M (k) = E m (k x ) + E n (k y ) + E q (k z ), (4) 

where E n (k) is the one-dimensional Bloch electron dispersion in the n-th band. The dipole 
matrix element is equal to [2] 



J fnq (x)X0mk ( X ) dx = ~ S m J IT S (<1 ~ k ) + X mn 8 {<i ~ k ) , where 

— OC J- 



2m r * , „ d 



X mn = — J U„ k (x) — U mk (x)dx. 
Q cell 



(5) 
(6) 



dk 



II. Crystal in Electric Field. Transfer Matrix Approach 

Consider an electric field S applied along the x-axis. The Schrodinger equation for motion of 

the particles with charge (-e) in the x-direction assumes the form: 



h 2 d 2 

H<p(x) = E<p(x) , where H = + U(x) + e£x . 

2m dx 



(7) 



where m is the bare electron mass, and potential U(x) = -a^S(x-na). Wave function (p(x) 

n 

in the interval (n - \)a < x < na between the crystal planes is that of a free electron in electric 
field S 



^,) = 4,Ai J-f,-A 



+ B Bi 



V A o 






(8) 



where A ft = 



( # v 3 



y2m eS j 



is the confinement length in electric field for bare electrons. 



Transfer matrix M for the transfer of the vector of coefficients through the period of the lattice 
in potential U(x) is defined as 



( A ^ 



= M. 



\ B nJ 



M=l + ^ 



Ai(Z„)Bi(ZJ Bi 2 (Z„) 

-Ai 2 (ZJ -Ai(Z„)Bi(Z„) 



(9) 
(10) 



where 1 is the unity matrix, / = 



2m a 



is an effective length of the periodic potential and 



z - m jtr-h. 



(in 



At the left boundary of the system n = n L = — N/2 + 1, and at the right boundary n = n R = 
N/2, where N = L/2a is the number of the unit cells corresponding to the thickness L of the 
nanofilm. The total transfer matrix M between the boundary planes of the nanofilm is given by a 
product of the known partial transfer matrices M n as 



n=NI2 

m= n M n 

n=-NI2+\ 



(12) 



This transfer matrix can, in fact, be calculated completely analytically using Mathematica 
software (by Wolfram Research). 

The coefficients of the wave function (8) at the left and right boundaries (denoted by the 
corresponding indices A and B) are related (transferred) as 



(A \ (A \ 
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= M 
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(13) 



Assuming the electrons are confined inside the nanofilm by high enough potential walls, we 
impose zero boundary conditions cp f ± -J = 0, which translate into 
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= 0. 
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V A oV2 



eS. 



= 0. 



(14) 



(15) 



These equations (13)-(15) form a homogeneous linear system of four equations for four 
quantities A , B , A , and B . The condition of existence of its nontrivial solution is the 

"L n L n R "R 

vanishing of its determinant, which, in principle, allows one to find the energy spectrum. Under 
this condition, one of these quantities can be set arbitrarily. We set, e.g., 



Aif-Lf-*-*^ 

A =1, B - ^ 2 eg -< , (,6) 



"I "I , . 

Bi 



V A ol 2 e ^ 



thus satisfying Eq. (14). The normalization of the wave function at the end of the procedure will 
produce the standard wave function. The remaining Eqs. (13) and (15) depend on energy E as a 
parameter. They can easily be resolved, resulting in a transcendental equation that we have to 
solve numerically. Note that except for this step, the rest of the problem is solved analytically, 
which dramatically increases numerical precision obtained at the end. Thus, in this case it is 
possible to reduce the Schrodinger equation, which is a partial-differential eigenproblem, to an 
explicitly analytical, though complicated, transcendental equation. As customary, we call this 
solution exact. (It is presently conventional to call analytical or exact a solution of a differential 
(eigen)problem that is reduced to explicit integrals, series, or a algebraic/transcendental equation. 
The entire solution, including their numerical evaluation, is often called semi-analytical.) 
In realty, to calculate the energy spectrum of the system, we have employed the shooting 
method, which can be described as the following. We set the trial value of E (this initial value 
should be negative enough). At the current step, we evaluate the transfer matrix of Eq. (12) and 
then coefficients A and B from Eq. (13). At the end of the current iteration, we compute the 

n R n R 

right-hand side of Eq. (15). Increasing the value of E, we find a point when this right-hand side 
changes sign, which correspond to the current eigenvalue of E. Then this procedure is repeated 
until required number of the eigenstates are evaluated. (This procedure reminds finding the right 
distance to target by gradually increasing the shooting angle, which gave it the name "shooting 
method".) With the spectrum of the system known, the wave functions are found analytically 
from Eq. (8). 



III. Crystal in Electric Field. Crystal Momentum Representation. 

A solution of Eq. (7) in the crystal momentum representation can be written as an expansion over 

the Bloch functions: 



JL : Li 

(p{x) = Y J \dka m (k)<?> mk (x). (17) 



m -nta 



Substituting Eq. (17) into Eq. (7) 

7i I a 

X jdka m (k)(E m (k)-E + eSx)<t> mk (x) = 

m -nla 

and after multiplying by <f nq {x) and integrating over x we arrive at [3] 

7T I a co 

(E„ (q) ~ E)a n (q) + eS ^ \ dka m ( k ) J C ( x ) x <P mk ( x ) dx = 



(18) 



m -Jtla 



(19) 



E n (q)~E + ieS — a n (q) + eSj^ a m (q)X mn = 0. 

Uq J m 



In one-band approximation (i.e. X mn = , if m ^ n ) Eq. (19) has the form 



E n (q)-E + ieS 



The solution to Eq. (20) is [3] 



d_ 

dq 



a n (q) = 0, where E = E - eSX ni 



(20) 



a n (q) = AQxpl — \(E n (q')-E)dq'>, where q is a constant, 



eS 



'/ (1 
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<P n ( X ) = A j" dk ex P1 ik x-— + — JE n (q)dq \u nk (x) 



eS 



. k 



(21) 



(22) 



IV. Carriers in Region of Stopping Points 

Consider eigenfunctions in the region of a stopping point x « (E - E n )/eS at the n-th band 



extrema at k = ±q , where E n = E n (±q ) . 



JL I u | 

q> n (x) = A [ dk Qxp< ik 



E~E n 
eS 



+ -^l(E n (q)-E n )dqiu nk (x) 



(23) 



In the spirit of the stationary phase approximation the biggest contribution to the integral of Eq. 
(23) comes from the interval of A: where {E n (q) - E n ) I eSa « 1. Therefore, the Bloch electron 

dispersion under the integral can be expanded as E n (q)» E n ±h 2 (q±q ) 2 /2m* , where m is 

the effective mass at k = ±q (with "+" in front of the quadratic term in the dispersion 

corresponding to electrons and "-" to holes). As a result for x « (E - E n )/e£ 

<P n { x )*u m (xy¥ n (E,x), (24) 



x ¥ n (E,x) = 2AJ dkcoslk 



E-E B 
eS 



h 2 (k-g ) 3 | 1 
2m eS 3 e£ 



\(E n (q)-E n )dq\ 



where ^(i?,*) is the envelope function. For the band extrema at the T-point q - the integral 
in Eq. (24) can be evaluated as 
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(25) 
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where A = 



' £ 2 



.1/3 



is the confinement length for a dressed carrier (electron or hole), and 



2m eS j 

C = const . Thus, we have shown that the envelope function in the region of a stopping point at 
the T-point is given by the Airy function 



( 1 f 



x ¥ n (E,x) = NAi 
where N is a normalization constant. 



+ 
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(26) 



V. Quantum Bouncing of Bloch Electrons 

Consider carriers that are confined between the walls of the film and r -point reflection position. 

Boundary condition at the wall X — —L/2 for conduction band quasi-electron envelope is 



W c (-L / 2) = , which results in condition 
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A~ 



L Ei-E 



e v 



+ 
2 eS 



he 



- d , where Ci is i-th root of 



Ai'(-x) . The energy of the i-th state at the conduction band edge can be expressed as 

Et=E cb -eS(L/2-CAe)- 



(27) 



The quasi-electron wave function can be written as <p ni (x) « Nu nQ (x)- Ai 



r \( L^ 



A 



x + - 



<i 



The 



v iv v ^y j 



normalization constant can be found following the approach of the Ref. [4] 



-LI2 



■( 1 
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r t\ 



x + 



L 



C, dx = A£ l Ai(-£ i ) 2 



(28) 



v ^ J j 



Finally, as follows from Eqs. (26), (27) and (28) the normalized wave-function of the quasi- 
electron has the form 



Pciix) 



1 



Ar^rMrc) 
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eg 



u c0 (x). 



(29) 



J) 



The intraband dipole matrix elements between the quantum bouncer quasi-electron states given 
by Eq. (29) can be also found following the recurrent approach of the Ref [4] 



(i \x\ j) = 8 H \-- + -A e C I + (1 - S„ ) - 



■A. 



(30) 



Boundary condition for quasi-hole envelope at the wall x-L/2 is given by v F n (Z/2) = 0, 



which leads to condition 



1 
A 



L Et-E v 
2 eS 



- E, t , where ^, is i-th root of Ai(-x) . The energy of 



the states is 

Et=E vt +eS{LI2-Z i \ h ). 

Normalization constant can be found from the integral [4] 



1/2 / 1 
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ft dx = AAi'{-$ i ) 2 
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As a result the normalized wave function of quasi-holes has the following form 



<P vi (x) 
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The dipole matrix element of the hole state in the quantum bouncer regime can be found as 



(Mj) =s « 



rh^h-^A 



v 



From Eq. (27) the energy of the conduction band edge can be found as 

E h (g) = E h -eg(L/2-£ 1 A e (S)). 

The valence band edge, as follows from Eq. (31), is given by 

E l (S) = E l +eg(L/2-^A h (g)). 

The band gap E (g) is thus equal to 

E g (S) = E g -eS(L-C l A e (S)-&Ak(S))- 



(34) 



(35) 



(36) 



(37) 



The metallization field g m found from condition E (£ m ) = is 






. (38) 

Since L » A e (g m ),A h (g m ) one can very well approximate expressions (35), (36), (37) and 
(38) by 

E ht (£) = E ht + e$L/2,E g (£) = E g - e£L and £ m =E g l eL (39) 



VI. Wannier-Stark regime 

If the state is confined between two reflections at stopping points it is disconnected from the 

walls. In this case the state is bulk-like [5]. Its normalized wave function is given by 



<Pn(x) 



/ sl/2 tela 



In J 



eg 



J dkexplik x-— +—\E n (q)dq\u nk (x) 



(40) 



The state is invariant with respect to transformation x — >x + d, E^>E + e£d . From the 

condition that the Bloch components at the Brillouin zone boundaries are equivalent, one obtains 
[3] the energies of the Wannier-Stark states 



(-E/eg) — + — \E n (q)dq = 2nl 



-7i i a 
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E = — E n (q)dq - eSal 



2k 



(41) 



(42) 
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We have found that the intraband dipole matrix elements between Wannier-Stark states are given 
by 



7i I a 



(I \x\ I') = — f dk E n (q) ex V {ika(l - /')} , for / * /' . 

9 77K> € J 
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-n I a 



(43) 



VII. Oscillator Strength and Dielectric Permittivity 

Dielectric permittivity of the nanofilm is given by the following expression 



s(co) = 1 + — y 



4x^(" f ~ n i )\( d x)if\ , An y-,Kn f -n i )\(d x ) if \ 



hv 



COy—CO — z'O 



1 + - 



fiLa' 



■Z- 



a>y -co- iO 



(44) 



Here n =Sn denotes population of a three-dimensional state, while n i is population of a state 
in x-direction. The factor 8 originates from 2 occupied bands in y- and z- directions; another 2 is 
due to spin degeneracy. Also (d x ) if =e(f\x\i) is a transition dipole moment corresponding to 

the transition frequency m if . The imaginary part of the dielectric function is 

T . . 32;r 2 ^ / .,.,. ,2., 

lme(co) = —^lS n i ~ n fW*W 5 ^ a if ~ ®) • ( 45 ) 

rl-LtCl if 

The combined oscillator strength in the low-frequency part of the spectrum is 

f ma ^ I i2 

-=«2>MI • (46) 

where the sum is taken over the low-frequency transitions with transition frequencies C0j and 

x j ={fW)- 

The imaginary part of dielectric permittivity averaged over the low-frequency band Aa> is 

T _, x 32tt e ^-ii, x |2 

lms{co) = J^L?^ x) '\ (47) 



./ 



The static dielectric constant is given by 

TiLa j cOj 

VIII. Additional Discussion: Breakdown, Possible Experimental Observations, and 
Some Technological Ramifications 

In this Section, we briefly discuss the possibilities for experimental observations and 
technological ramifications in more detail than in the Letter text. 

The metallization effect considered in the present Letter critically depends on the possibility to 
apply a rather high electric field S and observe the effect with or without damage (breakdown) 
to the crystal. The dielectric breakdown is a common name for a group of phenomena. The 
fundamental upper limit to the electric field that an insulator can withstand is imposed by the 
Zener tunneling, as described in the text of the Letter. There is also a number of phenomena 
called summarily as breakdown, which are due to acceleration of shallow-level electrons in the 
applied electric fields. Such shallow levels are produced by impurities (dopants), 



crystallo graphic defects, or may be surface states. Their acceleration by the applied electric field 
causes avalanche by impact ionization, resulting in the lattice damage. 

The breakdown fields are maximum in high quality insulators crystals. In particular, the high 
quality insulator nanofilms used as the gate insulators in field-effect transistors have thickness 
L~\ nm and work without a breakdown under potential difference Acp ^ IV corresponds to the 
internal fields S S; 0.1 V/A, which are just below the fields typical for the metallization 
threshold. It is actually somewhat likely that the breakdown field S BD ~ 0. 1 V/A in the gate 
insulators is defined by the present metallization effect. 

As we will discuss below in this Section, in the optically-applied fields, the metallization causes 
the population of the conducting states and, consequently, high values of the dielectric 
permittivity of the film, making it effectively a metal. The dielectric screening by the conduction 
electrons lead to much lower fields inside than outside. This is responsible for the high fields 
tolerated by the metals at their surfaces, cf. Ref [6]. The metallization causes a plethora of the 
well-known nanoplasmonic effects. In the strong fields, they can be observed with femtosecond 
resolution by various approaches of ultrafast nonlinear optics, including generation of harmonics 
[6], photoelectron emission [7-8], pump-probe spectroscopy, etc. In this case, even if the 
metalizing fields are strong enough to cause the breakdown and damage, the ultrafast (femto- 
and attosecond) measurements, in particular, the detection of harmonics or high-energy 
photoelectrons can be accomplished before the lattice disintegrates. 

The theory in the present Letter is developed in terms of the electric field S that is the internal 
(inside the material) electric field. Correspondingly, the experimental manifestations and the 
possibilities to observe this effect crucially depend on the way this field is created. While there 
may exist many approaches to do so, we will discuss two basic ways: (i) by the electromagnetic 
wave incident from infinity (the scattering formulation) and (ii) by means of two conducting 
electrodes that form capacitor plates between which the nanolayer of the material is placed. 

(i) This approach is typical for an experiment where an electromagnetic (optical, THz, 
microwave) wave is impingent on the nanolayer from a distant source in such a way 
that its polarization is normal to the nanolayer. In this case the field S of the wave 
inside the system is determined by the excitation field and the charges at the surfaces 
of the system induced by this field (the depolarization phenomenon). If the frequency 
of the field is far from any resonances of the system (which is true for a dielectric 



10 



nanofilm or a metalized nanofilm) then the field at the surface is equal to the field B x 
at the infinity (i.e., at the source) and the internal field is by a factor of £ less, where e 
is the dielectric permittivity of the nanolayer, which should be computed including 
the metallization effect itself. 

When the threshold field S m of the metallization is reached, the dielectric permittivity 
increases more than tenfold to s ^ 100, which leads to S = S^Js « S x . This means 
a strong dielectric screening, which tends to keep the internal field at values just 
above S m . Such a metal-like dielectric screening protects the crystal from damage 
and allows for the external field to be quite large. For instance, in experiments [6] the 
source field intensity was 7 W ~ 10 n W/cm 2 , while the system was resonant and the 
field at the surface had local intensity 1 ~ 10 12 -10 13 W/cm 2 . Such an intensity was 
not damaging for the metal of the nanostructure. Note that local intensity 



1 ~ 10 11 W/cm 2 corresponds to field S = J 0.1 V/A, which is a typical field 

for the metallization problem under consideration. 

It may also be of interest to apply much larger fields, which will be damaging. 
However, for films with thickness L~3 nm, as shown in the Letter, the required 
excitation pulse is only a few femtoseconds. In such a case, the metallization will 
occur and can be probed by similar ultrashort pulses before the nanosystem is 
destroyed. Such a type of experiments is characteristic of optics of ultrafast and 
ultrastrong fields. 

When expressed in terms of the internal field B , the metallization problem is a 
purely quantum-mechanical problem, which is considered in the Letter. However, 
expressed in terms of the source field S^ it is a deeply nonlinear problem due to the 
dramatic change of £ in the course of metallization. We will consider such a nonlinear 
problem elsewhere. 

Small (nanoscopic) defects of the nanofilm surfaces in the case of the metallization 
will lead to the appearance of nanoplasmonic hot spots just as it occurs at the nano- 
inhomogeneities of the metal surfaces - cf. Ref. [9]. While they may cause ionization 
and breakdown, they are of interest as a phenomenon. In contrast, the variations of 
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the thickness of the nanofilm under these excitation conditions will complicate effect 
but not lead to the breakdown. This is because the normal electric field in the 
nanofilm is S cx> /s irrespectively from the film thickness. Of course, sharp steps at the 
surfaces may create enhanced local fields as the defects do [9]. 

As we point out in the Letter, the large negative value of Re s and significant 
dielectric losses described by Im s are typical for plasmonic metals. The very 
existence the fields nanolocalized at the metal surfaces requires Re e < , which 
condition is satisfied for metalized nanofilms in a wide spectral area - see Fig. 3 (b)- 
(d). An important parameter in plasmonics is the so-called quality factor Q = 

6 Vim £ - ^ or me metauze d nanofilms, depending on the frequency and S this 
quality factor ranges Q~l — 10, which is about the same as for gold in the visible 
spectral range, 
(ii) The formulation of the metallization problem as an insulator crystal nanofilm 
between two conducting electrodes is that the defined internal field S is defined 
determined by ration of the external potential difference AU and thickness L. In such 
a case, the metallization, when developed may or may not lead to the breakdown 
depending on internal resistance of the source. If this resistance is low (in engineering 
terms, the source is a generator of potential difference), then the metallization will 
lead to very high currents and a damage of the crystal. Such a situation is typical for 
one of the most important potential applications of this theory: the supercapacitors. In 
such a case, the metallization is a defining factor for their energy capacity, which 
fundamentally limits the maximum field in the supercapacitor dielectric to the values 
S ~ 0.1 V/A, which by an order of magnitude lower than the limit imposed by the 
Zener breakdown. 

Another important application of this formulation is the gate insulator in field-effect 
transistors. In this case, also the metallization will lead to the breakdown of this 
insulator and damage of the transistor. The metallization imposes here a fundamental 
limit on the gate field, which is a defining factor for the transistor performance. Note 
that actual fields in the transistor-gate insulators are on the same order as the 
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metallization field S m -0.1 V/A. It is quite possible that the actual transistor-gate 
breakdown field S BD ~ 0.1 V/A is defined by the metallization. 
Considering, the variations of the nanofilm thickness, the metallization is defined by 
the potential drop across the nanofilm and not by the field 2 per se. This is evident 
the first of Eqs. (3) in the Letter, which states that the metallization field S m is 

defined as eS m L = E g , which implies that at the metallization point the total potential 
drop across the nanofilm AU = eg m L is fixed and equal to the bandgap E g . 

Therefore, in this case the variation in thickness will not lead to the breakdown in the 
thinnest part in contrast to other mechanisms of breakdown including the Zener 
tunneling and field-induced avalanche. However, the non-uniformity in the nanofilm 
thickness will lead to different limitations on the rate with field may change with the 
respect to the adiabaticity. Thinner films are more adiabatic and therefore more 
tolerant with respect to the metallization mechanism. One need also mention that the 
transistor (MOS) technology produces highly uniform oxide nanofilms used as the 
gate insulators, which will be very suitable for the purposes of studying the 
metallization effect. 
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